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ABSTRACT 

Analysis of the absorption lines in the afterglow spectrum of the gamma-ray burst 
GRB 010222 indicates that its host galaxy (at a redshift of z= 1.476) is the strongest 
damped Lyman-a (DLA) system known, having a very low metallicity and modest 
dust content. This conclusion is based on the detection of the red wing of Lyman-a 
plus a comparison of the equivalent widths of ultraviolet Mg I, Mg II, and Fell lines 
with those in other DLAs. The column density of H I, deduced from a fit to the wing 
of Lyman-a, is (5 ± 2)xl0^^ cm“^. The ratio of the column densities of Zn and Cr 
lines suggests that the dust content in our line of sight through the galaxy is low. This 
could be due to either dust destruction by the ultraviolet emission of the afterglow or 
to an initial dust composition different to that of the diffuse interstellar material, or 
a combination of both. 
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1 INTRODUCTION 

Gamma Ray Bursts (GRBs), the most powerful cosmic ex¬ 
plosions, are very important as probes of the high-redshift 
universe. This is due to several factors: (i) GRBs occur 
at cosmological distances, the most distant so far being at 
z=4.5 (GRB 000131, Andersen et al. 2000), (ii) they can 
be extremely bright and therefore could be seen (though 
for a very short period of time) up to redshift 20 (e.g. Wi- 
jers et al. 1998; Lamb & Reichart 2000), (iii) the burst and 
its afterglow illuminate the material situated between the 
GRB and the observer, providing an opportunity to study 
galaxies that otherwise would be too faint to be observed 
spectroscopically. This is particularly relevant for the host 
galaxies of GRBs where, in analogy to quasar (QSO) ab¬ 
sorption line studies, the strengths of the absorption lines 
can be used to measure important properties, such as their 
gas content, dynamical structure, and metallicity. In addi¬ 
tion, since they are transient phenomena, they do not cause 

* Based on observations made with the 4.2-m William Herschel 
Telescope, operated on the island of La Palma by the Isaac New¬ 
ton Group in the Spanish Observatorio del Roque de los Mucha- 
chos of the Institute de Astroffsica de Canarias. 
f e-mail: isabel@science.uva.nl 


the large scale “proximity effect” of quasars, i.e. they do not 
modify their environment at large distances, although they 
are expected to do so in their immediate vicinity (see sec¬ 
tion |^, (iv) the currently most popular physical explanation 
of this phenomenon is the so-called fireball model (cf. Rees 
& Meszaros 1992) which is a relativistic explosion caused by 
either the collapse of a massive star (Woosley 1993; Paezyn- 
ski 1998) or the merging of two compact objects (Lattimer 
& Schramm 1974). In the first case the rate of GRBs in the 
universe should be directly linked to the instantaneous star 
formation rate (SFR). This indeed seems to be supported 
by the observations, since the host galaxies of GRB^ are 
undergoing strong star-formation activity (e.g. Fruchter et 
al. 2001; Vreeswijk et al. 2001). Therefore, GRBs offer the 
potential to measure the cosmic star-formation rate (SFR) 
as a function of redshift. 

In this paper we report on spectroscopic observations 
of the GRB 010222 afterglow taken ~23 hours after the 
burst. We present evidence that its host galaxy is a damped 


1- Actually, this is true for bursts whose duration is longer that 2 
seconds, since they are the only ones which have been localized 
up to now. 
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Lyman-a absorption (DLA) system with the strongest col¬ 
umn density of neutral hydrogen and the lowest metallicity 
currently known. This has strong implications for both the¬ 
ories of the nature of GRBs and their host galaxies, and 
for the studies of DLA absorption systems. In particular, 
it shows that DLA systems with N(Hl) > 10^^ cm“^ do 
exist, and that the way to detect them may be by looking 
at GRB afterglow spectra rather than QSO absorption line 
spectra. This is of extreme importance because DLAs are 
key components of the universe at high redshift due to the 
fact that, although relatively rare, they account for most of 
the neutral gas available for star formation. 

In section ^ we present the observations and describe 
the calibration procedures. In section ^ the reduced spectra 
are analysed, resulting in a detailed identification of the ab¬ 
sorption lines, and the determination of the metallicity and 
dust content of the gas producing them. Finally, in section ^ 
we determine the nature of the galaxy hosting GRB 010222 
and compare it to the host galaxies of other GRBs, as well 
as with other Mg II and damped Lya absorption systems. 
The discussion and conclusions are presented in section ^ 
and section P, respectively. 

1.1 GRB 010222 

GRB 010222 was discovered on February 22, 2001, at 7:23:30 
U.T. by the Wide Field Gamera (WFG) on board ReppoSAX 
(Piro 2001). This burst is one of the brightest bursts ever 
recorded, with a fluence of (9.3 ± 0.3)xl0~® in the en¬ 
ergy range from 40 to 700 keV. Soon after its detection, 
an optical transient was discovered with coordinates RA = 
14‘'52“12!55 and DEC =-f43°0l'06'.'2 (J2000), and magni¬ 
tudes V~18 and R = 18.4±0.1 (Henden 2001; McDowell et 
al. 2001; Henden & Vrba 2001). A series of optical spec¬ 
tra of the afterglow taken shortly after the alert allowed 
the determination of the redshift of the host galaxy (Gar- 
navich et al. 2001; Jha et al. 2001; Bloom et al. 2001; Cas¬ 
tro et ah, 2001), which is 2 = 1.476. For a cosmology with 
Ho = 65 kms“^Mpc”^, Q,m = 0.3, A = 0.7, the lumi¬ 
nosity distance is D_l =3.6x 10^® cm, and the isotropic en¬ 
ergy output is 7.8x10®® erg (in ’t Zand et al. 2001), which 
makes GRB 010222 the third most energetic of all GRBs 
with known redshift. 

With regard to the host galaxy, observations of 
GRB 010222 at mm- and sub-mm-wavelengths have led to 
the interpretation that it is vigorously forming stars at a 
rate of ~500 Mq yr“® (Frail et al. 2001), whilst in the op¬ 
tical and near-infrared it is a faint (V=26.0±0.1) compact 
(FWHM~0'.'15) galaxy (Fruchter et al. 2001). 


2 WHT SPECTRA OF THE OPTICAL 
AFTERGLOW 

We obtained optical spectra of the afterglow of GRB 010222 
on 2001 February 23, at 06:03:08.976 U.T. (i.e. 22.66 hours 
after the burst) with the double spectrograph ISIS at the 4.2- 
m William Herschel Telescope (WHT), of the Isaac Newton 
Group (ING) situated at “El Roque de los Muchachos” ob¬ 
servatory, La Palma, Spain. We used an EEV12 and a TeK 
CCD (charged coupled device) as detectors, together with 
the B300B and R158R gratings in the blue and red arm, 


respectively. The total exposure time was 50 minutes. The 
spectral resolution, as measured from the width of the arc 
lines, is 3.3 A in the blue arm, and 5.8 A in the red arm. 
The nominal wavelength coverage for the blue arm is from 
2728 to 6267 A. However, the camera optics in the blue arm 
suffer from vignetting, resulting in a reduced effective wave¬ 
length range (3042 to 5954 A). The red arm covers the range 
from 6140 to 9115 A. The log of the observations is listed in 
Table 

The data were debiased, divided by a normalized flat 
field and extracted following standard procedures in IRAfQ 
The spectra were then calibrated in wavelength using Ar-Cu 
arc spectra. A comparison between the observed and theoric 
wavelength of the Balmer series in the standard star gives 
an average error of 1.15 A in the blue spectrum, and 0.91 
A in the red spectrum. Sky emission lines were removed by 
subtraction of the average sky background spectrum next 
to the object. However, in the region redwards of ~7600 A 
the sky substruction is difficult and remnant sky emission 
lines are left. This does not affect the subsequent analysis as 
we have not used this part of the spectrum in our analysis. 
Finally, the flux calibration was performed using standard 
star BD-|-26°2606. The atmospheric extinction was applied 
using the mean extinction curve for La Palma. We present 
the final spectrum in Figs. and 

3 ANALYSIS OF THE ABSORPTION LINES 
3.1 Line identification 

We identified several absorption lines in the spectrum of the 
GRB 010222 afterglow, belonging to (at least) three differ¬ 
ent absorption systems at redshifts z=1.476, z=1.156 and 
z=0.927 (Jha et al. 2001; Masetti et al. 2001; Castro et al. 
2001). Actually, the GRB 010222 afterglow is one of the 
richest in terms of line detections; a compilation of absorp¬ 
tion lines detected in GRB afterglow spectra can be found 
in Vreeswijk et al. (2001). In this paper we will concentrate 
on the system at the highest redshift, which we will adopt 
as belonging to the host galaxy of GRB 010222. 

The host galaxy’s absorption lines and their identifica¬ 
tion are compiled in Table 2, together with the rest-frame 
equivalent width (EW) and some atomic parameters (the 
oscillator strength, /, and the ionization potential of the ion, 
IP, in eV). For comparison, we have listed the EW published 
by Jha et al. (2001) and Masetti et al. (2001). The lines that 
we could not identify (neither at redshift 1.476, nor at the 
other two redshifts) are listed in Table 3. 

It is important to determine whether the lines are satu¬ 
rated or not, which is difficult given the low resolution of our 
spectrum. Most of the absorption lines are expected to have 
a doppler broadening less than a few tens of kms“®, which 
is much less than the spectral resolution (~ 220 kms“®). 
Therefore, we have compared the observed ratio of the line 
intensities to the theoretical one, which is given by the ra¬ 
tio of their oscillator strengths. Furthermore, the absorption 

§ IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1. Log of the spectroscopic observations of the afterglow of 
GRB 010222. The data were taken on 2001 February 23, about 22.66 
hours after the burst. The total exposure time is 50 minutes. 
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Instrument 

CCD/Grism 

Spect. range 

Spect. resol. 



(A) 

(A / kms“^) 

ISIS Blue arm 

EEV12 / B300B 

3042 - 5954 

3.3 / 220 kms-l 

ISIS Red arm 

TeK / R158R 

6140 - 9115 

5.8 / 230 kms-l 


""As measured from the FWHM of the arc lines 





2500 3000 3500 



Wavelength (Angstrom) 


Figure 1. The optical spectrum of the afterglow of GRB 010222. For each panel, the observed wavelength is given at the bottom, 
and at the top, the wavelength at the host galaxy rest-frame (i.e. corrected by red-shift z=1.476). The spectrum includes many, mostly 
saturated, absorption lines which can be identified with common resonance lines at the redshift of the host galaxy {z = 1.476) and at 
two other absorption systems at redshifts {z = 1.156 and 0.927). The telluric absorption lines are marked with ©. At the redshift of the 
host galaxy, the center of the Lyman-a line should be located at ~3011 A, which is just outside the wavelength range covered by our 
spectrum. However, we can see its red wing extending up to 3400 A (observed wavelength). We do not detect the interstellar graphite 
feature at 2175 A (rest-frame); the origin of the “bump” in the continuum around 7600 A is unclear. 
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Figure 2. The optical spectrum of GRB 010222. We zoom in those areas which contain the lines listed in Table ^ and we marked them 
with a dotted line. The signal-to-noise per resolution element is ~13. 

lines with an equivalent width larger than ~0.5 A are most 
likely saturated (see for example Pettini et al. 2001). Accord¬ 
ingly, most of the absorption lines detected in the afterglow 
spectrum of GRB 010222 are actually saturated. 


3.2 Detection of Lyman-a at z=1.476 

The strongest line in the spectrum of GRB 010222’s after¬ 
glow is actually Lyman-a, of which we can see its red wing 
extending up to ~3400 A (observer’s frame) or -1370 A 
(rest-frame). The center of Lyman-a should be at 3010.8 A 
in the observer’s frame, which is just outside the wavelength 
region covered by our spectrum, close to the atmospheric 
and instrumental cut-off. We have checked that this feature 
is not due to instrumental effects, such as flat fielding or 
vignetting. For this purpose, we have compared the spec¬ 


trum of GRB 010222 afterglow to that of the standard star, 
which was observed right after the exposure of GRB 010222 
was finished. We can see in Fig. ^ that the feature that we 
identify as the red wing of Lyman-a is absent in the stan¬ 
dard star spectrum and that the vignetting starts at wave¬ 
lengths shorter than —3000 A (observer’s frame) or —1210 A 
(rest-frame). This is what we expect according to the tech¬ 
nical documentation published by the ING group (Garcia 
Lorenzo 2000). We have multiplied both spectra by an arbi¬ 
trary constant, so they can be plotted in the same diagram, 
and binned them so that their general shape is more clear. 

Gonsequently, we are confident that the feature ob¬ 
served at the blue end of the spectrum of GBR 010222’s 
afterglow is the red wing of the damped Lyman-a absorp¬ 
tion line of the host galaxy. 
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Wavelength (Angstroms) 
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Table 2. Absorption lines identified at 2 = 1.476. The first column lists the observed wavelength, the second the identification, 
the third the redshift, and the fourth the equivalent width in the rest-frame of the host galaxy, i.e. divided by (l-|-z). The error on 
the equivalent width is estimated by varying the location of the continuum. The error in the wavelength position is ~1.15 A for 
the blue spectrum and 0.91 A for the red spectrum. For comparison, we list in the fifth and sixth columns the equivalent width 
published by Jha et al. (2001) and by Masetti et al. (2001). The next columns are the oscillator strength of the transition, / (taken 
from Bergeson & Lawler 1993; Morton 1994, 2001; Savage &: Sembach 1996; Bergeson et al. 1996; Vernor 1996 and Schectman et 
al. 1998), and the ionization potential of the corresponding ion, IP, in eV. Finally, in the last column, we mention whether the line 
is as well detected in absorption line systems towards quasars (e.g. Jannuzi et al. 1998), and two supernovae, SN 1987A (Blades et 
al. 1988) and SN 1993J (Marggraf & Boer 2000). 


Obs. A 

Identification 

Redshift 

Eq. Width 

EW Jha 

EW Mas 

/ 

IP 

detected as well 

(A) 



(A) 

(±4.9 h) 

(±lday) 


eV 


3129.2 

Sill 1260.42 

1.483 

1.54 ±1.18 



1.007000 

8.151 

qso, 87A, 93J 

3134.4 

Sill* 1264.74 

1.478 

1.74 ±0.20 



0.903400 

8.151 

87A 

3161.9 

Cl 1277.24 

1.475 

0.98 ±0.07 



0.096650 

0 

87A 

3223.3 

01 1302.17 

1.475 

0.94 ±0.03 



0.048870 

0 

qso, 87A, 93J 

3304.7 

CII* 1334.53“ 

1.476 

2.34 ±0.03 



0.1278 

11.26 

qso, 87A 

3449.8 

Si IV 1393.76 

1.475 

0.67 ±0.04 



0.5280 

33.49 

qso, sn87A, 93J 

3472.3 

Si IV 1402.77 

1.475 

0.44 ±0.01 



0.2620 

33.49 

qso, 87A, 93J 

3779.7 

Sill 1526.71 

1.476 

1.65 ±0.17 

1.4 ±0.4 

1.2 ±0.3 

0.12700 

8.151 

qso, 87A, 93J 

3832.4 

CIV 1548.20 

1.475 

0.52 ±0.10 



0.190800 

47.89 

qso, 87A, 93J 

3839.5 

CIV 1550.78 

1.476 

0.71 ±0.14 

1.9 ±0.4f 

2.3 ±0.2 

0.09522 

47.89 

qso, 87A, 93J 

3982.3 

Fell 1608.45*’ 

1.476 

0.67 ±0.15 



0.0580 

7.87 

qso, 87A 

4136.9 

AlII 1670.79 

1.476 

1.49 ±0.09 

1.1 ±0.3 

1.1 ±0.3 

1.800 

5.986 

qso, 87A 

4477.4 

Sill 1808.01“ 

1.476 

0.56 ±0.02 

0.5 ±0.3 

0.7 ±0.3 

0.00218 

8.151 

qso, 87A, 93J 

4592.2 

Aim 1854.72'* 

1.476 

0.73 ±0.09 



0.5390 

18.83 

qso, 87A, 93J 

4613.3 

Aim 1862.79 

1.477 

0.72 ±0.07 



0.2680 

18.83 

qso, 87A, 93J 

5011.2 

Znll 2026.14 

1.473 

0.21 ±0.24 

1.0 ±0.3 

0.6 ±0.2 

0.489000 

9.394 

qso, 87A, 93J 

5017.0 

Mgl 2026.48 

1.476 

0.62 ±0.02 



0.1120 

0 

qso, 87A, 93J 

5091.4 

Cril 2056.25 

1.476 

0.22 ±0.07 



0.1050 

6.766 

qso, 87A 

5107.1 

Znll 2062.66 

1.476 

0.31 ±0.03 

0.7 ±0.2 

1.1 ±0.2 

0.256000 

9.394 

qso, 87A 


-l-Crll 2062.23 





0.0780 

6.766 

qso, 87A 

5116.4 

CrII 2066.16“ 

1.476 

0.38 ±0.33 



0.05150 

6.766 

87A 

6380.2 

Mnll 2576.88 

1.477 

1.21 ±0.40 

0.7 ±0.2 

0.5 ±0.2 

0.3508 

7.435 

qso, 87A 

6404.4 

Fell 2586.65 

1.476 

1.42 ±0.68 

1.5 ±0.2 

0.8 ±0.3 

0.069125 

7.87 

qso, 87A, 93J 

6421.4 

Mnll 2594.50 

1.476 

0.57 ±0.38 

0.7 ±0.2 

0.8 ±0.3 

0.2710 

7.435 

qso, 87A 

6437.7 

Fell 2600.17 

1.476 

2.16 ±0.70 

1.9 ±0.2 

1.9 ±0.2 

0.2390 

7.87 

qso, 87A, 93J 

6922.3 

Mg II 2796.35 

1.476 

3.33 ±1.27 

3.0 ±0.2 

2.8 ±0.2 

0.6123 

7.646 

qso, 87A, 93J 

6941.7 

Mg II 2803.53 

1.476 

3.82 ±1.65 

2.7 ±0.2 

3.4 ±0.3 

0.3054 

7.646 

qso, 87A, 93J 

7063.6 

Mgl 2852.96 

1.476 

1.32 ±0.89 

0.9 ±0.2 

1.0 ±0.3 

1.81000 

0 

qso, 87A 


“ The line of CII 1334.53 A is broad (155 kms“^ after correction for the instrumental profile). It is also observed as a broad feature 
in the spectrum of SN 1987A (Blades, Wheatley, Panagia et al. 1998). 

The line of Fell 1608.45 A is identified by Masetti et al. (2001) as Znll 2062 A at 2 = 0.926, but in this case we should observe 
Znll 2026 A as well, which we do not. 

^ The Sill 1808.01 A line should be accompanied by another one at 1816.93 A, but the latter is contaminated by a sky emission 
line. 

Note that this line is identified by Masetti et al. (2001) as being Fell 2383 A at 2 =0.926. 

® This line is too strong which makes its identification uncertain. Maybe this line is blended with some other line. 

^ The separate doublet components are not resolved, hence the equivalent width corresponds to the sum of both components. 


3.3 Measuring the column densities 

The average redshift, based on the metal lines, of the 
absorption-line system with the highest redshift is 1.467 (we 
have not taken into account the blends of Znll+Cril). Ac¬ 
cording to Castro et al. (2001) there is some substructure in 
the gas forming these lines, with a separation of 106 kms“^. 
The spectral resolution of our spectrum does not enable us 
to resolve this, and therefore we can only deduce the global 
properties of the galaxy or region where the GRB has taken 
place. 

To measure the equivalent width, we used the routine 
splot within iraf. The line was integrated with respect to 
the local continuum, by just summing the flux over a given 


wavelength interval. No gaussian fitting was done, except 
in the few cases where the lines were partially blended, or 
the local continuum was noisy. In these cases, we checked 
the obtained EW by fitting a single Gaussian to the lines, 
yielding consistent results. To estimate the error on the EW 
measurement, we varied the position of the local continuum. 
The final value of the EW is the average of all measurements, 
while the error is conservatively taken as the difference be¬ 
tween the highest and the lowest values obtained. The re¬ 
sults are given in Table where we have listed the EW 
published by Jha et al. (2001) and Masetti et al. (2001) too. 
In most cases, the three independently obtained values are 
consistent within the errors. However, in two cases the val¬ 
ues do not agree: (a) the ClV doublet at 1548, 1550 A which 
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Table 3. Unidentified lines. The 
equivalent width is as measured in 
the observer’s frame. 


Observed A 

(A) 

Equivalent Width 

(A) 

3138.3 

2.69 ±0.18 

3178.5 

2.32 ±0.07 

3410.9 

2.27 ±0.31 

3709.8 

0.78 ±0.03 

4389.6 

0.57 ±0.17 

4488.6 

0.53 ±0.07 

4518.1 

0.58 ±0.03 

4525.2 

0.32 ±0.00 

4633.9 

1.49 ±0.26 

4681.2 

0.63 ±0.01 

4856.7 

0.91 ±0.46“ 

4859.6 

(sum) 

4914.0 

0.45 ±0.03 

5100.3 

0.55 ±0.02 

7340.1 

7.02 ±0.11 

7519.9 

1.59 ±0.27 

8504.5 

7.63 ±0.37 


“ This line is blended with the line 
at 4859.6 A. 


can be due to the fact that we resolve it, while Jha et al. 
(2001) and Masetti et al. (2001) do not; (b) the Znll + CrII 
2026 blend, but again this can be explained by a difference 
in the spectral resolution: if we add the equivalent width of 
Mg I 2026 A line to that of the Zn/Cr blend, then the values 
agree well. 


3.3.1 Hi column density 

We have deduced the column density of HI by fitting a Voigt 
profile to the red wing of Lyman-a. For that purpose, we cor¬ 
rected the spectrum for redshift, and normalized the contin¬ 
uum to unity. Subsequently, we produced model profiles of 
damped absorption lines, for three values of N(H I): 3x 10^^, 
5x10^^ and 7x10^^ cm“®, and compared them with the 
observed Lyman-a wing. Even though the observed spec¬ 
trum is noisy and only part of the line profile is covered, 
a large value of N(Hl) is required to match the observed 
red damping wing. The value we have adopted is N(Hl) = 
(5±2) xlO^^ cm“®. This is formally (i.e. within the errors) 
consistent with the value deduced from X-ray observations, 
(2.5±0.5)xl0^^ cm"^ (in’t Zand et al. 2001). 


3.3.2 Column density of the other lines 

As we have discussed earlier, most of the other lines seem 
to be saturated, and are, therefore, not suited to calculate 
their column density. We limit ourselves to those lines that 
are very likely not saturated. These are the Zn II and Cr II 
lines, given their intrinsic low cosmic abundance and the 
small value of the oscillator strength. Possibly the lines of 
Fe II 1608 A, and Mn II 2594 A are unsaturated too, but nev¬ 
ertheless, we regard the abundances determined from these 
transitions as lower limits. 


In order to deduce the equivalent widths of the indi¬ 
vidual lines in the Zn + Cr blend, we start by measuring 
the EW of Znll 2026 A which is mainly dominated by Zn^, 
and the EW of Cr II 2056 A which is not blended. The total 
EW of Znll + CrII 2026 A is 0.21, but the contribution 
of Cr II 2026 A can be neglected, because the ratio of the 
oscillator strengths of CrII 2056 A to CrII 2026 A is ~22. 
Therefore, we assume that the EW of Znll 2026 A is 0.21 A. 
This implies that the EW of Znll 2062 A is ~0.11 A. Since 
the total EW of the Znll -|- CrII 2062 A blend is 0.31 A, 
then the EW of the CrII 2062 A becomes ~0.2 A. The EW 
of the Cr II 2066 A line has a very large error (87% relative 
error). Therefore, this line has not been taken into account 
to deduce the column density of Cr'*'. 

Once we know the EW, we calculate N^oi using the fol¬ 
lowing relationship, valid in the limit of weak lines, between 
the (rest-frame) equivalent width and the column density 
(SchefHer & Elsasser 1987): 

12 EW 1 

iVcoi = 1.131 X 10^^ X -^-7 cm, 

where / is the oscillator strength of the transition, and A 
the wavelength in cm. 

In the subsequent analysis, we have taken for each ion 
(Zn II, Cr II, Fe II, and Mn II) the average Ncoi of the individ¬ 
ual transitions, if more than one unsaturated line could be 
measured. We assume that the abundance of each element 
X is closely approximated by X’'" because the second ionisa¬ 
tion potential is above 13.6 eV. The results are summarized 
in Table ^ and discussed in the following sections. 


4 THE HOST GALAXY OF GRB 010222: 
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Figure 3. Evidence that the host galaxy of GRB 010222 is a damped Lyman-a system. Left: Our blue spectrum clearly shows the red 
wing of Lyman-a, which is matched with three Voigt profiles corresponding to a column density of neutral hydrogen of 3x10^^, 5x10^^ 
and 7x10^^ cm~^. Before the fitting, the spectrum has been corrected for redshift and the adjacent continuum normalized to unity. 
Right: The spectrum of GRB 010222 (top) is shown next to that of the standard star (bottom); the feature that we identify as the red 
wing of Lyman-Q is not present in the spectrum of the standard star. Both spectra have been multiplied by an arbitrary constant so 
that they can be displayed in the same panel, and they have been binned so that their general trend is more obvious. 


A DAMPED LYMAN-a SYSTEM 

An absorption system is classified as a damped Lyman-a 
(DLA) system if the colunm density of neutral hydrogen is 
larger than 2x10^° cm“^. The column density of HI in the 
host galaxy of GRB 010222 is N(Hl) = (5±2) xlO^^ cm“^. 
This, by definition makes the host galaxy of GRB 010222 a 
damped Lyman-a system. Further support to the classifica¬ 
tion of the host galaxy as a DLA comes from the analysis of 
Mg I, Mg II, and Fell lines. According to Rao & Turnshek 
(2000), if the absorption system has EW(MgII 2796) > 0.6 
A and EW(FeII 2600) > 0.5, it has 50% probability of being 
a DLA system. This is obviously the case for GRB 010222, 
as already pointed out by Masetti et al. (2001). Further¬ 
more, all DLA systems in the Rao & Turnshek sample have 
Mg II 2796, 2803 A doublet ratios below 1.5, as is the case for 
GRB 010222 (Fig. |^. The relationship between EW(MgII 
2796) versus EW(FeII 2600) also indicates that the host 
galaxy of GRB 010222 is a DLA: On the one hand both 
measured EWs are larger than 0.6 A, and on the other hand 
their ratio is ~1.5, very close to the average value of 1.45 
deduced by Rao & Turnshek (2000) (Fig. |). 

Furthermore, the equivalent widths of other low ioniza¬ 
tion lines in the spectrum of GRB 010222 are quite large 
compared to those of absorption-line systems present in 
QSO spectra. For example, the EW of Sill 1526 A is 1.65 
±0.17 A, while the values measured in the sample of Steidel 
& Sargent (1992) range from 0.02 to 1.87 A. Equally strong 
are the lines of Znll 2026,2062 A and CrII 2026,2056,2062 
A. If we compare for example their EW with those of the 
data collected by Pettini and collaborators (1994, 1997 and 
1999), we see, again, that GRB 010222 has the highest val¬ 
ues (even though it has the lowest metallicity, see later). 

In order to estimate the metallicity and dust-to-gas ra¬ 
tio we compare the abundances of elements that are known 
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to be little depleted onto dust to that of refractory elements 
that are easily depleted on to dust grains. The best elements 
for such purpose are Zn (not depleted) and Cr, Mn, Ee and 
Ni (large depletion). Therefore, Zn is a good indicator of the 
metallicity, whilst Cr, Mn, Fe and Ni can be used to mea¬ 
sure the dust content (for a more extended analysis of why 
these elements are useful to study the metallicity and dust 
content, see for example Pettini et al. 1990). 

We adopt the usual notation, i.e. 


[X/H] 


Yeol(A) Yeol(X) 


log5(X),(l) 


where 5{X) is the depletion of element X, or in other words, 
the fraction of element X in gas form, relative to its solar 
meteoritic abundance. The latter values have been taken 
from Savage & Sanders (1996). 

For the zinc abundance, and thus a measure of the en¬ 
richment of the gas by metals, in the host of GRB 010222 
we find a value of [Zn/H] = -2.3 ±0.5, indicating that the 
host galaxy has a very low metallicity, about 200 times less 
than the solar fractional Zn abundance. If we adopt N(Hl) 
= 2.5x10^^ cm“^ (in’t Zand et al. 2001) then [Zn/H] be¬ 
comes -2.01, about a factor 100 below solar, thus the host 
galaxy would still be amongst the most metal poor DLAs 
known. The lowest metallicity claimed up to now in DLAs 
is [Zn/H] = -2.07 ±0.1 (Molaro et al. 2001). 

The dust content of the host galaxy can be deduced 
by comparing the Zn abundance to the abundance of some 
refractory elements: 


[X/Zn] = log{ 


S{X) 

S{Zn) 


log 


/gas ± fdu 


( 2 ) 


where X is Cr, Fe or Mn, S{X) its gas-phase abundance, and 
/gas and /dust are the fraction of element X in gas and dust 
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respectively. The value of [X/Zn] = -0.3 corresponds to the 
case where the fraction of X in dust is equal to the fraction 
of X contained in gas form (/gas = /dust). The results are 
summarized in Table The abundances of Zn, Cr and Mn 
are consistent with each other, while that of Fe is slightly 
smaller. 

A few words of caution before reaching any conclusion: 
The iron line Fell 1608 A has an EW slightly larger than 
0.5 A, and therefore could be saturated, so its abundance 
should be regarded as a lower limit. The same applies to 
Mn, which in addition has a nucleosynthetic pattern that 
is not well understood, so we may be seeing effects of nu¬ 
cleosynthesis, and not of depletion on dust (Pettini et al. 
2000). Also, it is important to realize that there might be 
effects due to the ionization of the gas caused by the local 
UV radiation field of OB stars and/or of the GRB itself. 
These effects would raise the value of [Zn/Cr] which could 
mimic the absence of dust (Howk & Sembach 1999). How¬ 
ever, further observational proof of the lack of dust in the 
host of GRB 010222 is the absence of the “2175 A dip”, 
which is usually attributed to graphite dust grains and de¬ 
tectable down to Av ~ 0.1 mag (see also Jha et al. 2001 and 
Masetti et al. 2001). 

In order to compare the host of GRB 010222 with DBAs 
detected in QSO spectra, we have compiled the data pub¬ 
lished by Pettini and collaborators (1994, 1997 and 1999). 
In Fig. M we plot the column density of neutral hydrogen, 
N(Hl), versus the Zn abundance [Zn/H] (i.e. the metallic- 
ity) and the abundance of Gr'*' (indicator of the dust con¬ 
tent) versus that of Zn'*'. Clearly, the host of GRB 010222 
is exceptional in a number of ways. First, the low value of 
the metallicity is what one would expect from the extrap¬ 
olation of the apparent correlation between N(Hl) and the 
metallicity. There is some controversy regarding the reality 
of such correlation. Dust bias in DLA selection is one of the 
possible explanations for the absence of high column den¬ 
sity, metal rich absorbers. A more plausible explanation is a 
cross-sectional bias against intercepting the center of galax¬ 
ies. For a more detailed discussion we refer the reader to El¬ 
lison et al. 2001 (and references therein). Second, the abun¬ 
dances of Zn and Cr are the lowest in the whole sample. The 
correlation of Zn and Cr abundances is caused by the fact 
that they are both Fe-peak elements and are produced by 
similar nucleosynthetic reaction pathways. In Galactic stars, 
Zn and Cr actually track each other very closely, hence any 
departure from [Cr/Zn]=0 indicates depletion onto dust. If 
there was no depletion on to dust we would expect a one- 
to-one ratio in Fig ^ Instead, we observe a scatter of points 
below this line which is due to the presence of dust. 


4.1 Other cases of DBAs in host galaxies of GRBs 

This is not the first time that a GRB host galaxy has 
been confirmed as a damped Byman-Q system. Jensen et 
al. (2001) detected Byman-a in the afterglow spectrum of 
GRB 000301C, and deduced a neutral hydrogen column 
density \ogN{HV) = 21.2±0.5 cm“^, qualifying the host 
galaxy of GRB 000301C as a DBA. Another case where the 
line of Byman-Q is as well detected is the host galaxy of 
GRB 000926, where the column density of neutral hydrogen 
is \ogN{H\) ~ 21.3 cm~^ (Fynbo et al. 2001). Although 


both cases have N(Hl) lower than in GRB 010222, they are 
still at the high end of the range found for typical DBAs. 

It is useful to compare these results with those of 
Galama & Wijers (2001). From the analysis of the X-ray 
spectra of GRBs afterglows, they find evidence for high col¬ 
umn densities (xi3xl0^^ cm”*^) of circumburst material. 
The X-ray analysis has the advantage that it measures the 
column density of all atoms, irrespective of whether they 
are in gas or dust. However, at 1-2 keV rest-frame energies, 
the X-rays only measure the abundance of elements whose 
K edge is at or below that energy (roughly elements with 
Z<10) and therefore one needs to know the metallicity in 
order to deduce N(HI). The usual practice is to assume so¬ 
lar metallicity which means that the real N(Hl) could be 
even larger, supporting the idea that GRBs always happen 
in DBA systems]^ 

Furthermore, analysis of the equivalent widths of Mg II 
2796,2803AA, Mg I 2853 A and Fell 2600 A lines provides 
further support to the idea that the sites of GRB explosions 
are DBA systems. We have taken from the literature the 
EWs of those lines seen in GRB afterglow spectra and added 
these points to the diagrams (figs. 24 and 25) presented in 
Rao & Turnshek (2000). The results are shown in Fig. ^ 
It is evident that not only the host of GRB 0102222 would 
qualify as a DBA on the basis of these diagrams, but also 
the hosts of GRB 970508, GRB 990123, GRB 000926, and 
possibly of GRB 990510 and GRB 991216 (Metzger et al. 
1997; Andersen et al. 1999; Kulkarni et al. 1999; Gastro et 
al. 2001; Vreeswijk et al. 2001). These results reinforce the 
idea that all GRB host galaxies are DBAs. In addition, Rao 
& Turnshek (2000) conclude, apart from the fact that all the 
systems in their sample with EW(MgI 2582) > 0.7 A are 
DBAs, that there is no particular correlation between the 
EWs of Mg II 2796 A and Mg I 2582 A lines. However, this 
changes if we add the values for the GRB mentioned before: 
there is a clear trend for both EW being proportional to 
each other (see Fig. 


5 DISCUSSION 

From the analysis of the optical spectrum of the afterglow 
of GRB 010222, we conclude that its host galaxy is a DBA 
absorption line system, with the strongest column density 
of neutral hydrogen and the lowest metallicity ever known, 
and with low dust content. On the other hand, based on 
the bright millimetre and sub-millimetre constant emission 
(3.74 ±0.53 mjy at 350 GHz), Frail et al. (2001) conclude 
that the host galaxy of GRB 010222 is a dusty galaxy with 
an intense burst of star formation (~500 MQyr”*^) while 
in the optical and near-infrared it is a blue sub-luminous 
galaxy. 

There is an apparent contradiction, because Pettini et 
al. (1999), based on the lack of metallicity evolution from 
z=0.5 to z=3.5, concluded that DBAs do not trace the bulk 

^ In principle, when calculating N(HI), one should account for 
any intervening galaxies. Typically, those values will be similar to 
our own galaxy at high latitude, a few times 10^** cm“^, which 
means that if one finds such high values, N(HI) > lO^*^ cm”*^, it 
is most likely due to the host galaxy. 
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Figure 4. The host galaxy of GRB 010222 has the highest column density of HI and the lowest metallicity of the whole sample of 
DLAs. The data of the latter have been taken from Pettini et al. (1994, 1997 and 1999). The values of the abundances of Zn and Cr 
and of the column density of HI are what one would expect from an extrapolation of the correlations seen in both diagrams. A word 
of caution however: the correlation between N(HI) and [Zn/H] is thought to be due to selection effects; there is no physical reason why 
there should not be absorption systems with a large N(HI) and high metallicity (see more details in text). The dotted line in the right 
diagram corresponds to the one-to-one correlation between the abundances of Zn and Cr that we would observe in the total absence of 
dust. Filled squares correspond to DLA systems with low dust content and open squares to systems with high dust content. The filled 
pentagon represents the host of GRB 010222. The errors in the Pettini sample vary from 1 to 20 per cent for N(HI), and from 0.2 to 0.4 
dex for the abundances of Zn and Cr. 


Table 4. Abundances and dust content in the host galaxy of GRB 010222. The column densities of Zn and Cr are obtained by averaging 
the results for the individual transitions Znll 2026,2062 AA and CrII 2056,2062 AA respectively. The column density of Mn and Fe 
has been obtained from a single transition, Mn II 2594 A and Fe II 1608 A. We assume that the abundance of each element can be 
approximated by the abundance of its first ionised ion, since the second ionisation potential is in all cases greater than 13.6 eV. For 
comparison, we list as well the range of abundances of these elements through different lines of sight in our Galaxy (see Savage & Sembach 
1996 and references therein). Solar meteoritic abundances are taken from Savage & Sembach (1996). 


Ion 

Ncol 

[X/H] 

[X/Zn] 

fdust 

fgas 

[X/H] 

[X/H] 

[X/H] 

[X/H] 


cm“^ 

010222 

010222 

Halo 

Disk+Halo 

Warm Disk 

Cool Disk 

Zn+ 

(1.25±1.64)xl0l3 

-2.3±0.4 

_ 

- 


rsjQ 



Cr+ 

(6.21±2.46)x1013 

-2.6±0.7 

-0.34 

~1.2 

(-0.38,-0.63) 

(-0.72,-0.88) 

(-1.04,-1.15) 

(-2.08,-2.28) 

Mn+ 

>(3.53±2.35)x1013 

>-2.7±0.9 

-0.43 

~1.7 

^0.47,-0.72) 

^0.66) 

^0.85,-0.99) 

(-1.32,-1.45) 

Fe+ 

>(5.05±1.13)xl0i'‘ 

>-3.5±0.6 

-1.26 

~17 

(-0.58,-0.69) 

(-0.80,-1.04) 

(-1.19,-1.24) 

(-2.09,-2.27) 


of star-forming galaxies at these redshifts. Indeed, low red- 
shift imaging studies (e.g. Le Brun et al. 1997) have shown 
that DLAs represent a wide range of morphological types. 
At higher redshifts, the low spin temperatures inferred from 
the 21 cm absorption line studies indicate that DLAs rarely 
represent massive galaxies and may be more analogous to 
metal-poor dwarfs (Chengalur & Kanekar 2000). On the 
other hand, the properties (column density, profile shapes, 
velocity dispersion and number of components) of Mg II lines 
seen in the spectra of QSOs seen at z ~1.5 are very simi¬ 
lar to those observed in r] Carina nebula. This means that 
these lines could be associated with an expansion event, like 
superbubbles driven by star formation (Danks 1999). How¬ 
ever, Bond et al. (2001), after comparing the kinematics and 
redshift evolution of strong Mg II absorbers, conclude that 


superwinds and star-forming galaxies can account for a sub¬ 
stantial fraction of these, but not for DLAs which show dif¬ 
ferent kinematics and no red-shift evolution. 

Furthermore, DLA systems having N(Hl) > 10^^ cm“^ 
are very rare. In fact, to our knowledge, the host of 
GRB 010222 is the only one known. One of the explana¬ 
tions often suggested was that such systems would have a 
large dust content and therefore would be missed in opti¬ 
cal surveys, but this cannot be the dominant effect because 
the survey of radio selected QSOs (and therefore circum¬ 
vent the dust issue) by Ellison et al. (2001) did not find 
any DLA with a large N(Hl). Another possible reason pro¬ 
posed by Schaye (2001a) is that clouds having such high 
column densities have a large fraction of molecular hydro¬ 
gen, and such clouds are much smaller and shorter lived 
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(adapted from Rao & Turnshek, 2000) 




Figure 5. Are all GRB host galaxies damped Lyman-u systems? Left: The ratio of the equivalent width of the two Mg II doublet lines 
is below 1.5 for all the confirmed DLA systems in the sample of Rao &: Turnshek (2000) and for those GRB for which this ratio has 
been measured (see label in figure). Right: Another indication that GRB host galaxies are DLAs: The equivalent width of MgII 2796 A 
is correlated to that of Fell 2600 A for the DLA systems of the sample of Rao Sz Turnshek (2000). The dotted line corresponds to 
the correlation given in Rao Sz Turnshek (2000). The GRB host galaxies for which this information is available follow this correlation. 
Moreover, all GRB hosts have equivalent widths of Fell 2600 A larger than 0.8 A, and of MgII 2796 A larger than 2 A. When the 
individual components of the Mg II lines were not be resolved, two points are drawn in the diagram connected by a vertical line: the 
upper value is the total measured EW and the lower value is the EW divided by two. 



Figure 6. Relationship between the equivalent width of Mg II 2796 A and that of Mg I 2852 A. There is a clear trend for both quantities 
being correlated, that is not evident if one takes only into account the DLA systems of the Rao Sz Turnshek (2000) sample. The dots 
joined with a vertical line correspond to those cases where the Mg II doublet lines were blended: the upper value is the total measured 
EW and the lower value is the EW divided by two. 


than clouds with low molecular fractions. Furthermore, at 
high redshift, there is a large population of compact galax¬ 
ies with high star formation rates and large-scale outflows 
(called Lyman Break Galaxies, LBG), whose spectra resem¬ 
ble those of DLA systems seen through the line of sight to¬ 


wards QSOs. These sight-lines will most likely intersect the 
outer parts of the galaxy, while the LBGs will probe the in¬ 
ner parts of the galaxy (for a detailed discussion see Schaye 
2001b). The cross-section of such clouds will be very small, 
therefore making the chance bisection of a background QSO 
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unlikely. If GRB 010222, and GRBs in general, are associ¬ 
ated with young, violent star formation regions, that is, with 
giant molecular clouds, which most likely have large N(Hl), 
then the place to look for high column density DBAs is the 
spectra of the GRB afterglows, instead of those of QSOs. 

Another issue of controversy is dust. To explain the 
fact that the optical/UV emission of the host galaxy is faint 
compared to the sub-milimetre emission. Frail et al. (2001) 
propose obscuration by dust. However, if the system had 
a high content of dust, there should be a higher amount 
of extinction and a larger depletion of certain elements, in 
particular Cr, than the one observed (see section One 
possible explanation is that the dust is destroyed by the UV 
and X-ray radiation from the GRB and its afterglow. Ac¬ 
cording to Fruchter et al. (2000), the X-ray emission of the 
burst can destroy dust (by grain charging mainly) within 
the first minutes out to 20 pc. Also, according to Waxman 
& Draine (2000), the UV flash of the burst can sublimate 
dust in the first hours after explosion. This would explain 
the low extinction and the absence of the 2175 A feature 
which is due to small graphite particles (see also Galama 
& Wijers 2001 and references therein). However, according 
to Draine (2000) if the GRBs are associated with a molecu¬ 
lar cloud, the UV flash would vibrationally excite H 2 which 
would then cause an absorption at wavelengths < 1650 A 
which would be visible over tens of days. We do not observe 
such feature in our spectrum. Moreover, if dust is destroyed 
by the GRB and its afterglow emission, the ionization front 
should have affected the relative abundances of the different 
ions for the same element. Another factor that we must take 
into account is that the initial size distribution of dust parti¬ 
cles could be different to that in the diffuse Galactic media, 
with a dehciency of small grains and an increase of large 
ones. This is what is expected in very dense enviromnents 
-such as molecular clouds- where small particles are system¬ 
atically removed (Kim 1994) and it is indeed observed in the 
Small Magellanic Cloud (Weingartner & Draine 2001, Welty 
2001) and in some Seyfert galaxies (Maiolino et al. 2001). If 
this applies to the environment where the GRBs take place 
then the feature due to silicate grains at 9.7 /rm should be 
absent as well. 

A word of caution: the low value of [Zn/H] that we 
hud for GRB 010222 does not necessarily mean that its 
host galaxy is especially metal poor. It could just mean 
that cold, dense, giant molecular clouds are particularly de¬ 
pleted. Given the fact that the QSOs are unrelated to the 
absorption-line sites and the GRBs may not be, it is unclear 
whether the GRB host galaxies are special or only the GRB 
environments within hosts. The low resolution of our spectra 
does not allow us to disentangle these two cases. There is one 
case, GRB 000926, where Castro et al. (2001) find that the 
circumburst material is formed by two different clouds with 
a velocity separation of 168 kms“^. In these two clouds the 
relative abundances of Zn and Cr are similar to each other, 
thus indicating that the metallicity and dust depletion are 
characteristic of the host galaxy as a whole, and ruling out 
the hypothesis of grain destruction in the vicinity of the 
burst. In the case of GRB 010222 there is also some sub¬ 
structure with a velocity separation of 106 kms~^(Castro 
et al. 2001). If these subsystems had similar abundances, 
then the metallicity and dust content that we Hnd would 
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be representative of the galaxy as a whole, and again, dust 
destruction by the burst itself would be ruled out. 


6 CONCLUSIONS 

The observational evidence points out that GRB 010222 has 
taken place in a damped Lyman-a system having the highest 
column density of neutral hydrogen and the lowest metallic¬ 
ity known up to now, and with low dust content. This result 
reinforces the idea that a small cross-section, and not dust, 
is the cause of not finding DLA systems having very large 
N(Hl). If the system associated with GRB 010222 is not 
atypical, this means that we may be missing an important 
fraction of HI at cosmological distances which in turn may 
have an impact on Hdla (i.e. the fraction of neutral gas 
in DBAs as a function of the closure density), although the 
exact contribution will depend on their cross section, which 
is unarguably small. 

It is not clear whether the dust is destroyed by the GRB 
itself, or if its size distribution is different, being devoided of 
small particles. In the latter case the feature due to silicate 
grains at 9.7 /rm should be absent also. 

Finally, there is strong evidence that other GRB host 
galaxies (or at least their locations) are as well damped 
Byman-a systems. Since they seem as well associated with 
star formation, we conclude that they may occur in a sub¬ 
class of DBAs that is associated with active star formation. 
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